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Abstract 

, A massive binary black hole (BBH) is inevitably formed in a merged galactic nucleus before the black 

' holes finally merge by emitting the gravitational radiation. However, it is still unknown how the BBH 

I evolves after its semi-major axis reached to the sub-parsec/parsec scale where the dynamical friction with 

_ ■ the neighboring stars is no longer effective (the so-called the final parsec problem). In this paper, we 

Qj^', propose a new mechanism by which the massive BBH can naturally coalesce within a Hubble time. We 

i study the evolution of the BBH with triple disks which are composed of an accretion disk around each black 

hole and one circumbinary disk surrounding them. While the circumbinary disk removes the orbital angular 
momentum of the BBH via the binary-disk resonant interaction, the mass transfer from the circumbinary 
disk to each black hole adds some fraction of its angular momentum to the orbital angular momentum 
of the BBH. We find that there is a critical value of the mass-transfer rate where the extraction of the 
' orbital angular momentum from the BBH is balanced with the addition of the orbital angular momentum 

to the BBH. The semi-major axis of the BBH decays with time whereas the orbital eccentricity of the BBH 
^ I grows with time, if the mass transfer rate is smaller than the critical one, and vice versa. Its evolutionary 

^ I timescale is characterized by the product of the viscous timescale of the circumbinary disk and the ratio 

c/3 , of the total black hole mass to the mass of the circumbinary disk. Since a minimum value of the critical 

^ ' mass-transfer rate is larger than the Eddington accretion rate of massive black holes with masses in the 

IO^Mq to 10^ Mq range as far as the evolutionary timescale is shorter than a Hubble time, it is promising 
' that the critical mass-transfer rate is larger than the mass transfer rate. Most of massive BBHs, therefore, 

^ I enable to merge within a Hubble time by the proposed mechanism, which helps to solve the final parsec 

problem. 
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1. Introduction the common gravitational potential due to the dynami- 

cal friction with neighboring stars. When the separation 
. Massive black holes in galactic nuclei are con- between two black holes becomes less than 1 pc or so, an 
sidered to have co-evolved with their host galaxies angular momentum loss by the dynamical friction slows 
(Ferrarese&Merritt 2000; Gebhardt et al. 2000; Magorrian down due to the depiction of the stars on orbits inter- 
^ ■ et al. 1998). Since galaxies are well-known to evolve secting the BBH (Saslaw et al. 1974). This is the second 
JH I through frequent mergers, this strongly suggests that evolutionary stage. Finally, the BBH coalesces rapidly 
black hole growth is mainly caused by black hole mergers if the semi-major axis decreases to the point where the 
and subsequent accretion of gas (Yu & Tremaine 2002; Di emission of the gravitational radiation becomes an effi- 
Matteo et al. 2005). If so, a massive binary black hole cient mechanism to remove the orbital angular momen- 
(BBH) is inevitably formed before the black holes merge tum of the BBH. The transition from the second stage 
by emitting the gravitational radiation. Even if there are to the final stage is, however, considered to be the bottle- 
transiently triple massive black holes in a galactic nucleus, neck of evolutionary path for the BBH to coalesce because 
the system finally settles down to the formation of the of cutting off the supply of the stars on intersecting or- 
massive BBH by merging of two black hole or by eject- bits. This is called the final parsec problem (see Merritt 
ing one black hole from the system via a gravitational & Milosavljevic 2005 for a review). 

slingshot (Iwasawa et al. 2006). Recent hydrodynamic Many authors have tackled the final parsec problem in 

simulations showed the rapid BBH formation within sev- the context of the interaction between the black holes 

eral Gyrs by the interaction between the black holes and and the stars, but there has been still extensive discus- 

the surrounding stars and gas in gas-rich galaxy merger sions (Roos 1981; Makino 1997; Quinlan & Hernquist 

(Mayer et al. 2007). 1997; Milosavljevic & Merritt 2001; Makino & Funato 

It is widely accepted that massive BBH mainly evolves 2004; Merritt et al. 2007; Sesana et al. 2007). There is 

via three stages (Begelman et al. 1980; Yu 2002). Firstly, other possible way to extract the energy and the angular 

each black hole sinks independently towards the center of momentum from the BBH by the interaction between the 
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Fig. 1. An artist's impression of a massive BBH with triple 
disks on a parsec/subparsec scale of merged galactic nucleus. 
The BBH is surrounded by a circumbinary disk (CBD), from 
which the gas transfers to the central binary, and then the 
accretion disk is formed around each black hole. 

black holes and the gas surrounding them. In some cases 
of circular binaries with extreme mass ratio, the secondary 
black hole could be embedded in the gas disk around the 
primary black hole and then migrates to the primary black 
hole. This kind of the binary-disk interaction could also be 
the candidate to resolve the final parsec problem (Ivanov 
et al. 1999; Gould & Rix 2000; Armitage & Natarajan 
2002; Armitage & Natarajan 2005). 

Hayasaki et al. (2007) found that if the BBH is sur- 
rounded by the gaseous disk with the nearly Keplerian 
rotation (i.e. the circumbinary disk: CBD), the gas can 
be transferred from the CBD to each black hole. The 
mass transfer leads to the formation of the accretion disk 
around each black hole (Hayasaki et al. 2008), and then 
the BBH system finally has the triple disks which are com- 
posed of the accretion disk around each black hole and 
one circumbinary disk surrounding them (see Fig. 1 for a 
schematic view of the BBH system). There is, however, 
little known how the BBH evolves in such a triple-disk 
system. We, therefore, study the evolution of the mas- 
sive BBH interacting with the triple disks. Although we 
mainly discuss the case of a massive BBH system, our 
results can be applied to other possible context such as 
the compact binaries, the young binary star formation, 
and the extrasolar planet formation because of the scal- 
ing nature. The plan of this paper is organized as follows. 
In Section 2, we describe the derivation of basic equations 
governing the orbital evolution of the BBH and the formu- 
lation of the gravitational torque of the BBH which acts on 
the CBD. The solutions for their evolutionary equations 
are, then, reported in Section 3. Section 4 is devoted to 
discussions, and finally we summarize results in Section 5. 

2. Basic equations 



third law, by which the BBH is defined. The total energy 
of the BBH E'b is written by 



2a 



(1) 



where a is the semi-major axis of the BBH and /.j = 
il/iM2/Mbh is the reduced mass of the BBH. Here, il/bh = 
M1 + M2 is total black hole mass. Mi is the primary black 
hole mass, and M2 is the secondary black hole mass. By 
differentiating both sides of above equation, the energy 
dissipation rate of the orbital motion can be obtained as 
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(2) 



where Mi is the mass accretion rate of the primary black 
hole, M2 is the mass accretion rate of the secondary black 
hole. The dot over each physical quantity shows the time 
differentiation, unless otherwise noted. 

The orbital angular momentum of the BBH Jb is writ- 
ten by 



Jb = ^a^r^b VT~~e^, 



(3) 



where e is the orbital eccentricity, Ob is the angular fre- 
quency of the BBH. The relationship between the orbital 
energy Ey, and the orbital angular momentum Jb can be 
written, from equation (1) and (3), 

Ob Jb 
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The change rate of the orbital angular momentum Jb can 
be written as 
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where Mbh = Mi + M2 is the total mass-accretion rate. 
The orbital evolution of the BBH is generally governed 
by equations (l)-(5). Since the total mass accretion rate 
is much lower than the black hole mass, the second term 
and the third term in the right hand-side of equations 
(2) and (5) can be neglected. Actually, Hayasaki et al. 
(2008) showed by the Smoothed Particle Hydrodynamics 
simulations that the mass accretion rate is two orders of 
magnitude lower than the mass transfer rate. Therefore, a 
set of the basic equations of the orbital evolution is finally 
given by 



a 
a 

and 



Eb 
El, 



l-e2 
respectively. 
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2.1. Interaction between the disk and the binary black 
hole 



We assume that two black hole are gravitationally 
bounded and their motions are followed by the Kepler's 



We consider the resonant interaction between the BBH 
and the CBD. The CBD is assumed to be geometrically 
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thin, be aligned with the orbital plane of the BBH, and 
have nearly Keplerian rotation with no self-gravitation. 
The binary potential can be expanded by Fourier double 
series, 



(8) 



where m is azimuthal number. I is the time-harmonic num- 
ber. The potential component 0m./, can be written by 

(i(f7bi) / <^cos{m9 -m],t), (9) 
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which give rise to a corotational resonance at radius where 
ilp = f2 and the Lindblad resonances at radii f2p = ± 
{n/m). Here, Qp is the pattern frequency of the binary 
potential which is defined by 



ftp = — fib, 
m 



(10) 



At is the epicyclic frequency, and the upper and lower signs 
correspond to the outer Lindblad resonance (OLR) and 
the inner Lindblad resonance (ILR), respectively. The 
radii of these resonances for the CBD are given by 



''CR = -r 
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and 
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Here, the epicyclic frequency k ~ for the nearly 
Keplerian disks. The standard formula for toques at the 
LR is given by Goldreich & Tremaine (1979) 



LH,_m(m±l)^^S(A~2m)^0^, 
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(13) 



where A = In (/),„/ /d In r|LRs- If m > 1, A = — (to + 1) for 
the OLR in the CBD (cf. Artymowicz & Lubow (1994)). 
Similarly, the torque at the CR is written as 



^CR^2m3^!Ml 



(14) 



Note that the angular momentum is added from the bi- 
nary to the CBD via the LRs and the CR. On the 
other hand, the viscous torque formula derived by Lin 
& Papaloizou (1986) is written as 



Tvis = 37rQ!ssS17^r'* 



H 



(15) 



where ags is the Shakura-Sunyaev viscosity parameter 
(Shakura & Sunyaev 1973). Since the viscous torque re- 
moves the angular momentum from the CBD, the CBD is 
truncated and forms the gap between the CBD and the bi- 
nary if the viscous torque is less than the resonant toque. 
The inner edge of the CBD is formed at the given reso- 
nance radius where the viscous torque is balanced with 
the resonant torque: 
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where the summation is taken over all combination of (m, 
1) which give the same resonance radius. Actually, cri- 
terion (16) is determined only by the viscous torque and 
the torques from the OLR of the lowest-order potential 
component, because the torques from the OLR dominate 
those from the CR in the CBD and high-order potential 
components contribute little to the total torque, even if 
the eccentricity of the binary is not small (Goldreich & 
Tremaine 1979; Artymowicz et al. 1991; Artymowicz & 
Lubow 1994). Following Artymowicz & Lubow (1994), 
the inner edge of the CBD rjn is approximately determined 
by {m,l) = (2, 1) resonant torque in the binary with a low 



eccentricity: rjn = {m+ l/l) 



2/3, 



2.08a. It is confirmed 



by Hayasaki ct al. (2007) that this value is also applica- 
ble for the BBH with the moderate eccentricity e = 0.5. 
Therefore, the BBH transfers most angular momentum to 
the CBD via the 1:3 resonant radius nn ~ 2.08a. Below, 
the set of {m,l) is regarded as (2,1), unless otherwise 
noted. 

3. Orbital evolution of massive binary black holes 

The long-term evolution of the semi-major axis and of 
the orbital eccentricity are mainly driven by the resonant 
interaction between the BBH and the CBD. In this sec- 
tion, we firstly describe the evolutionary relation between 
the semi-major axis and the orbital eccentricity, the evo- 
lution of the semi-major axis, the evolution of the orbital 
eccentricity, and finally the effect of the mass transfer from 
the CBD on the orbital evolution of the BBH. 

3.1. Relation between the semi-major axis and the orbital 
eccentricity 

As the motion of the BBH slowly vary with time by res- 
onantly interacting with the CBD, there is an adiabatic 
invariant which can be regards as the orbital angular mo- 
mentum of the BBH (Landau & Lifshitz 1960). In such 
the system, the energy dissipation rate of the BBH is pro- 
portional to the product of the change rate of the adiabatic 
invariant and the characteristic frequency of the system, 
e.g., rip(cf. Lubow & Artymowicz 1996): 

Eh = flpX. (17) 

From equations (4), (10), and (17), equation (6) is rewrit- 
ten as 



(18) 



a m Jb 

which determines the evolution of the semi-major axis of 
the BBH. Using above equation, equation (7) is rewritten 
as 



l-e2 



l-^yi 



(19) 



From equation (18) and (19), we can obtain the differential 
equation relating a to e: 
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Fig. 2. Evolution of the semi-major axis a as a function of 
the orbital eccentricity e. The semi-major axis, a, is normal- 
ized by the initial value of the semi-major axis ao- The initial 
value of the orbital eccentricity is indicated by eo . The solid 
line, the dashed line, and the dotted line are the a-e relations 
of eo = 0.0, of eo = 0.5, and of eo = 0.9, respectively. 
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a m yj\ 

The above equation can be integrated, and then, a can be 
expressed as the function of e: 
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ao 



1 



1 



1 



(21) 



where oq and eo are the initial value of the semi-major axis 
and the orbital eccentricity, respectively. Figure 2 exhibits 
the evolution of the semi-major axis a normalized by ao 
as a function of the orbital eccentricity e. It is noted from 
the figure that the semi-major axis decays with time as 
the orbital eccentricity grows with time, and vice versa. 
When the BBH evolves towards a black hole merger, the 
eccentricity growth is faster than the orbital decay because 
the orbital eccentricity is more rapidly close to 1.0 than 
the semi-major axis is close to 1.0. 

3.2. Evolution of the semi-major axis 

Assuming that the CBD completely absorbs the orbital 
angular momentum of the BBH and is quasi-steady state, 
the change rate of the orbital angular momentum Jb can 
be written as Jb — ~ Jcbd, where Jcbd is the change rate 
of the angular momentum of the CBD. 



= -2 



^ Jcbd 
m Jb 



(22) 



Since the orbital angular momentum is always transferred 
from the BBH to the CBD via the inner edge of the CBD 
Tin, Jcbd can be approximately estimated as 
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where q is the mass ratio of the secondary black hole to 
the primary black hole, the scope of which is 0.1 ^ q < 
1.0 because the CBD is not truncated and thus the triple 
disk model is broken down if q is less than 0.1. Also, 
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S, and r^?^ is the viscous timcscalc of the 
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(24) 



where the initial value of a semi-major axis ao = Ipc, 
and the CBD is assumed to be isothermal for simplicity. 
Equation (22) is finally written as 

(25) 

a tc\aoJ 

Here tc shows a characteristic timescale of the evolution 
of the BBH, which is defined by 



to 
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where r^f (if* is the viscous timescale of the CBD when 
ao, and C is independent of both the semi-major axis 
a and time t: 
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Figure 3 represents the characteristic timescale tc of the 
massive BBH evolution as function of the black hole mass 
in unit of the solar mass. Here, we set the Tq = IQ'^K 
and ass = 0.1, which are typical values for a disk in active 
galactic nuclei (AGNs). The ratio of the CBD mass to 
the black hole mass is set as McBD/^bh = 10^^, which 
ensures that the CBD is stable for the self-gravitation be- 
cause Toomre's Q value is Q = (McBD/Mbh)"^(Jf/r) > 1 
(Toomre 1964). The solid line denotes tc in the range 
Mhh = 10^-10^ A/0. As seen in Fig. 3, tc is the longer as 
the black hole is more massive. 
Integrating equation (25), 

-=(l-f)', (28) 
ao V tc J 

which determines the evolution of the semi-major axis of 
BBH. Figure 4a shows the evolution of the semi-major axis 
a normalized by the initial value of the semi-major axis ao. 
It is noted from the figure that the semi-major axis rapidly 
decays with time regardless of the orbital eccentricity, and 
the BBH finally merges at t — tc- 

3.3. Evolution of the orbital eccentricity 

Substituting equation (25) into equation (20) can be 
written as 
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Fig. 4. (a) Evolution of a semi- major axis normalized by the initial value of the semi-major axis, ag, of the BBH for < t < tc- 
The solid lino shows the evolution of the semi-major axis. (6) Evolution of an orbital eccentricity of the BBH for <t <tc. The 
solid thin line, the dashed line and the dotted line show the evolution of the orbital eccentricity in the case of an initial value of the 
orbital eccentricity eo = 0.0, eo = 0.5 and eo = 0.9, respectively. 
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Fig. 3. Characteristic timescale, tc, of massive BBH evolu- 
tion as a function of the total black hole mass in units of the 
solar mass, M^i^/Mq, in the case of AfcBD/A^bh = 10~^, the 
Sakura-Sunyaev viscosity parameter ogg =0.1, the black-hole 
mass ratio M2/M1 = 1.0, and the disk temperature Tq = 10* A". 
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which determines the evolution of the orbital eccentric- 
ity. Integrating both sides, we can obtain the evolutionary 
timescale of the orbital eccentricity 
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The evolution of the orbital eccentricity e is 



e=v/(l + ^W)(l-^W), 
where rj is 



vm = - 



m 
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(30) 



(31) 



(32) 



Figure 4b represents the evolution of the orbital eccentric- 
ity of the BBH in the case of cq = 0.0, e = 0.5, and e = 0.9, 
respectively. This figure clearly shows that the orbital ec- 
centricity grows with time and is close to 1.0 within t = tc 
regardless of the initial value of the orbital eccentricity. 
The higher value of the initial eccentricity is more rapidly 
close to 1.0. 

3.4^. Effect of the mass transfer from the circumbinary 
disk 

Hayasaki et al. (2007) found that the mass transfer from 
the CBD to each black hole occurs every binary orbit. 
Since the transfered mass has the relative angular mo- 
mentum to each black hole, the accretion disk is formed 
around each black hole (Hayasaki et al. 2008). In this sec- 
tion, we investigate the effect of the mass transfer on the 
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Fig. 5. (a) Same format as Fig. 4a. The solid line, the dashed line, and the dotted line show the evolution of semi-major axis when 
(Mt) /^crit = 1-25, {Mt)/ Merit = 2.0, and (MT)/A'4rit = 0.5, respectively, (ft) Same format as Fig. 46. The soUd thin Une, the solid 
thick line, the dashed line, and the dotted line show the evolution of the orbital eccentricity when (eo, (Mx)/A^crit) = (0.9,1.25), 
(eo,{MT> /Merit) = (0.9,2.0), (eo, (MT>/M„it) = (0.5,1.25), and (eo, (MT)/M,rit) = (0.0,0.5), respectively 



evolution of the massive BBH. 

The averaged torque, which is added to the accretion 
disks by the mass transfer, during one orbital period is 
given by 

'/^l + g)2(A/T) Jb 



1 



(33) 



where (Mt) denotes the averaged mass-transfer rate dur- 
ing one orbital period, which can be expressed as the sum 
of the averaged mass-transfer rate to the primary black 
hole (the primary transfer rate, (Afx.i)) and the averaged 
mass-transfer rate to the secondary black hole (the sec- 
ondary transfer rate, {Mt,2))- (Mt) = (Mt,i) + (Mt,2>. 
Unless otherwise noted, we call (Mr) the mass transfer 
rate. The accumulated mass around each black hole ac- 
cretes onto each black hole by which the angular momen- 
tum is transferred outward and is finally added to the 
binary due to the tidal torque (cf. Kato et al. 2008). The 
torque added to the binary by this process, Jadd, can be 
estimated as 



pbhd " 



bhd7 



(34) 



where Jpbhd and Jsbhd are the torque added from the ac- 
cretion disk around the primary black hole to the BBH 
(the primary torque) and the torque added from the ac- 
cretion disk around the secondary black hole to the BBH 
(the secondary torque), respectively. By defining the ratio 
of the primary transfer rate to the secondary transfer rate, 
qx ~ (.^^t.2)/(Mt,i), the primary torque can be written 



by 
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^ ms, 1 
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where Ph is the orbital period of the BBH, rvis,i is the 
viscous timescale measured at the circularization radius 
rc,i of the accretion disk around the primary black hole, 
and rins,i is the radius of the marginally stable circular 
orbit of the primary black hole. Here, we approximate 
t/pbhd by that r^s,! is much smaller than Tc^i. Similarly, 
the secondary torque can be written by 



sbhd ■ 



qriMT) Ph 



1 



QT Tvis,2 



^/GM^, 



(36) 



where rvis,2 is the viscous timescale measured at the cir- 
cularization radius rc,2 of the accretion disk around the 
secondary black hole. Substituting equations (33) and 
(36)-(35) into equation (34), we can obtain 



Jadd = JJt- 

Here, / is the parameter which is defined as 



(37) 
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(38) 



where fb is the orbital velocity of the BBH, Cs,i and Cs.2 
are the sound velocity of the primary disk and the sound 
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velocity of the secondary disk, respectively. For simplic- 
ity, we assume that Cs,i~Cs,2 and the ratio of Cs.i to Vh is 
constant during the evolution of the BBH. The parameter 
/ shows how much the torque of the mass transfer is con- 
verted to the torque of the BBH. When / equals to 1.0, the 
torque of the mass transfer is completely converted to the 
torque of the BBH. The value of / is, however, substan- 
tially much smaller than 1.0 because the orbital velocity is 
much faster than the sound velocity at the circularization 
radius in general. 

The angular momentum balance in the BBH system can 
be expressed as 



Jh — — ^CBD + ^add 

A (Mt) / a 

Jb 




-1/2 



(39) 



and thus, from equation (25), the differential equation of 
the semi-major axis can be written as 



tc 




-1/2 



(40) 



where Mc^a is the critical mass transfer rate which can be 
defined as 

3McBD 



Merit = 



f_CBD 
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Here, we use equation (26) and (27). Integrating both 
sides of equation (40), wc obtain 
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Figure 5a displays the evolution of the semi-major axis 
of a massive BBH with triple disks. It is noted from the 
figure that the semi-major axis a decays with time when 
Merit > (Mt), while the semi-major axis increases with 
time when Merit < (Mt) • In the hmit of (Mt) /Merit < 1 -0, 
the dotted line corresponds to the solid line of Fig. 4a. The 
growth rate of the semi-major axis is more steep when 
(Mt)/ Merit = 2.0 than when (Mt) /Merit = 1-25. This 
means that the semi-major axis more rapidly grows with 
time as accretion disks are more massive. 

Similarly, the differential equation of the orbital eccen- 
tricity e is modified by the mass transfer as follows: 
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By integrating both sides of above equation, wc can obtain 
the evolutionary equation of the orbital eccentricity. 



e=V(l + ^W)(l-'/W), 
where rj is 
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Figure 5& represents the evolution of the orbital eccen- 
tricity of a massive BBH with triple disks. The fig- 
ure shows that the orbital eccentricity decays with time 
when {Mt) / Merit > 1-0, whereas it grows with time when 
{Mt)/ Merit < 1-0. This is independent of the initial 
value of the orbital eccentricity eg. The rate of evolu- 
tion is more rapid as {Mt)/ Merit is larger than 1.0 when 
{Mt)/ Merit > 1 and as {Mt)/ Merit is smaller than 1.0 
when {Mt)/ Merit < 1. In the limit of {Mt)/ Merit < 1-0, 
the dotted line corresponds to the solid line of Fig. 46. 

Which conditions arc more promising, (MT)/Merit < 1 
or {Mt)/ Merit > 1 ? The answer is {Mt)/ Merit < 1. In 
order to confirm this answer, we consider the minimum 
value of the critical mass transfer rate in what follows. If 
te is longer than a Hubble time, the binaries never merge 
within a Hubble time. We, therefore, confine our argu- 
ment to the case in which te is shorter than a Hubble 
time. Since the possible value of te is, then, a Hubble 
time at a maximum, the minimum value of the critical 
mass transfer rate can be estimated as 

-1/3 



Me 



m ( m + 1 



1 M 



bh 



in 



(l + g)2 I \ I J /, 

where in = I/Hq ~ 1.37 x 10^°yr is a Hubble time, f, 
is the maximum value of / which can be expressed as 



(46) 



: 27rass 



1 



Here Wb,min is 



GMbh 



-1/3 



7I/2 



Cs,l 
^b.mii 



2(1 + g) 



(47) 



(48) 



where ah = G^/Aa^ is the hardening radius where the 
binding energy of the BBH exceeds the kinetic energy of 
the stars, and a is the ID velocity dispersion of the stars in 
the galactic nucleus (Quinlan 1996). Assuming a approx- 
imately equals to the observed stcller velocity dispersion, 
we obtain from the best current estimate of M-tr relation 
(Merritt & Milosavljevic 2005), 
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with /3= 1/(4.86 ±0.43). Substituting equations (47)-(49) 
into (46), we can obtain 

2.0 X 10^1 /m\ / g \ / 0.1 \ 

j) lloJ 

2/3+1 r; 



Me; 



1.37(1.66 ±0.24)2% 
/200kms-i" ^ 



Cs,l 



(io^Mq 



Mq 



,(50) 



On the other hand, the Eddington accretion rate A/Edd is 



Ah 



47rGA/bh 



Edd 



0.2 



(. 



A/, 



bh 



Ma 



(51) 



ckt VIO^A^q/ [ yr 

where c is the velocity of light and kt is the Thompson 
opacity coefficient. The ratio of Ajfciit.min to A/Edd can be 
written as 

/r.,^,^^ _1 \ 2 



Ah 



Edd 



/200kms~i 
1.37(1.66 ± 0.24)2% 1^ 



A/, 



(lO^AfQ 



2^; 



(52) 



rmin 
^crit 



— ) (-)• 

ass/ Vl.O/ 

lOkms"^ in the outer region 
is always larger than A^Edd 



Since the sound velocity Cg^i 
of a typical AGN disk, 

for supermassive black hole with the mass in the IQ'^Mq to 
lO^A^Q range. Hayasaki et al. (2007) actually confirmed 
that the mass transfer rate is sub-Eddington rate if the 
gas is supplied to the CBD at the Eddington rate. It 
is, therefore, promising that most massive BBH evolves 
towards the orbital decay with the growth of the orbital 
eccentricity. 

4. Discussion 

We study how a massive BBH with triple disks evolves 
through the binary-disk interaction. In the long-term, the 
BBH resonantly interact with the CBD at the OLR ra- 
dius via which the orbital angular momentum is mostly 
extracted by the CBD. In the short-term, the gas is trans- 
ferred from the CBD to each component of the BBH. This 
process is repeated every binary orbit during the evolu- 
tion of the BBH, by which the accretion disk is formed 
around each black hole, and then each accretion disk vis- 
cously evolves. Some fraction of the angular momentum 
of each accretion disk is finally converted to the orbital 
angular momentum of the BBH due to the tidal torque 
of each black hole. When the mass transfer rate corre- 
sponds to a critical value, no BBH evolves because the 
angular momentum exchange caused by these binary-disk 
interactions is balanced. 

The direction of BBH evolution is determined by the 
critical mass transfer rate. If the mass transfer rate is 
smaller than the critical mass-transfer rate, the separation 
between two black holes decays with time whereas the or- 
bital eccentricity grows with time and is finally close to 
1.0. The grown-up eccentricity will, then, cause the black 
hole to plunge into the CBD, resulting in the enlarge- 
ment of the gap between the central binary and the CBD 
(cf. Artymowicz & Lubow 1994). This direct interaction 
would make the semi-major axis shorten because of the 



friction between the black hole and the gas in the CBD. 

What happen if the mass transfer rate is large enough 
to be over the Eddington accretion rate? Since the mass 
transfer rate is, then, larger than the critical mass trans- 
fer rate, the semi-major axis rapidly increases with time 
whereas the binary orbit become circularized with time. 
When the binary separation reaches to the hardening ra- 
dius, the growth of the semi-major axis will be again 
stalled by the dynamical friction with the neighboring 
stars. Thus, the BBH is long-lived in this case. 

Where is the orbital angular momentum of the BBH 
finally going? The angular momentum absorbed by the 
CBD is gradually transferred outward with the gas in the 
CBD. The gas around the BBH will, therefore, be dis- 
persed if no gas is supplied to the CBD. However, if the 
gas is supplied to the CBD during the BBH evolution, the 
outwardly transferred angular momentum will be removed 
by a supply source. This is essentially same problem as 
how the gas supplies to AGNs. Again, if the gas is sup- 
plied to the CBD at any rate, the CBD is likely to be 
a quasi-steady state. Such a state is, however, consid- 
ered to be different from the steady state of accretion disk 
around a single black hole. The orbital evolution of the 
BBH should, therefore, be solved with the evolution of 
the CBD. We will tackle this problem by the numerical 
simulation in the forthcoming paper where both the full 
torque of the BBH and the evolution of the CBD will be 
taken into account. 

Figure 6 shows the coalescent timescale of the mas- 
sive BBH with A/bh = IO^A^q and A/cBo/Afbh = 10"^. 
Two black holes firstly evolve towards the shrinkage of 
the binary separation due to the dynamical friction with 
the neighboring stars (Begelman et al. 1980). Next, the 
orbital evolution is stalled at Oh where the evolutionary 
timescale by the dynamical friction is over a Hubble time 
even if the emission of the gravitational wave dissipates 
the binding energy of the BBH. However, the BBH can 
track a new evolutionary path due to the binary-disk in- 
teraction. The new evolutionary path is denoted by the 
solid thick line of Fig. 6. In addition, the rapid growth 
of the orbital eccentricity allows the connection to the fi- 
nal evolutionary path where the BBH evolves towards the 
coalescence by emitting the gravitational radiation, since 
the timescale of the final evolutionary path in an eccentric 
binary is the factor (1 — e2) ''/2 times shorter than that in a 
circular binary (Peters 1964). The combination of the or- 
bital decay and the eccentricity growth makes it possible 
for the massive BBH to coalesce with a Hubble time. 

One of the most interesting problem is how many mas- 
sive BBHs are there in the Universe. It is noted from Fig. 3 
and Fig. 6 that the BBH more rapidly coalesces as the 
black hole mass is less massive and the CBD is more mas- 
sive. This suggests that a supermassive BBH with masses 
more than ~ 1O^^^A^0 is long-lived as a binary. Also, the 
relatively less massive BBH with masses lO^-'^A^Q would 
sequentially merge in early Universe and grow to a more 
massive, single BH. The host galaxies with such the grown 
black holes would suffer a major merger, resulting in a su- 
permassive BBH formation. If so, we can infer, by using 
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Fig. 6. Evolutionary timescalo of a massive BBH with triple disks from lOOpc scale to 10~^pc scale. The total black hole mass is 
A/bh = 10* Mq with equal mass ratio M2/M1 = 1.0, the ratio of the CBD mass to to the total black hole mass is Mcbd /A^bh = 10~^ , 
and the ratio of mass transfer rate to the critical mass transfer rate is {Mx) / Merit =0.5. The dashed line shows the first evolutionary 
track in which two black holes get close each other by their angular momentum loss due to the dynamical friction between the black 
holes and the stars surrounding them. The dotted line shows the third evolutionary track in which the BBH finally coalesce by the 
dissipation due to the emission of the gravitational radiation. The thick line shows a new, second evolutionary track in which the 
angular momentum of the BBH is removed by the binary-disk interaction. The horizontal dot-dashed line shows a Hubble time 
tu ~ 1.37 X lO^Oyr. 



M-a relation, that such supermassive BBHs are observed 
in well-developed galaxies with the gas-rich nuclei. 

5. Conclusions 

For the purpose of resolving the final parsec problem, 
we study the evolution of a massive BBH with triple disks 
on the parsec/subparsec scale. Our main conclusions are 
summarized as follows: 

(1) . The orbital evolution of the BBH is characterized 

(2) . There is a critical mass-transfer rate, which deter- 

mines the evolutionary direction of the BBH with 
triple disks. When the mass transfer rate is smaller 
than the critical one, a semi-major axis decays with 
time whereas the orbital eccentricity grows with 
time. The orbital eccentricity is rapidly close to 1.0 
during ic/2 even if the BBH is initially on a circular 
orbit, while the semi-major axis is 0.0 during tc re- 
gardless of the orbital eccentricity. The high value 
of the orbital eccentricity results in a rapid merger 
within a Hubble time through the emission of the 
gravitational radiation. 

(3) . When the mass transfer rate is larger than the criti- 

cal one, a semi-major axis grows with time whereas 
the orbital eccentricity decays with time. In such 
a system, the semi-major axis reaches to the hard- 
ening radius of the BBH, typically Oh ~ Ipc when 
Mbh = IO^Mq, and its orbit is completely circular- 
ized. 



(4). Since the minimum value of the critical mass- 
transfer rate is larger than the Eddington accretion 
rate of a massive black hole with mass in the IO^A/q 
to IO^A/q range as far as the evolutionary timcscale 
is shorter than a Hubble time, it is plausible that the 
mass transfer rate is smaller than the critical one. 
Most BBH, therefore, prefers to evolve towards a 
rapid coalescence within a Hubble time. 
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